Large Excitonic Effects in the Optical Properties of Monolayer M0S2 
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The band structure and absorption spectrum of monolayer M0S2 is calculated using the Go Wo 
approximation and the Bethe-Salpeter equation (BSE), respectively. We find that the so-called A 
and B peaks in the absorption spectrum arise from strongly bound excitons (0.7-0.8 eV) localized 
in distinct regions of the Brillouin zone and not from a split valence band as commonly assumed. 
Furthermore, we find the minimum band gap to be of the indirect type. This seems to conflict 
with recent experimental results showing stong luminescence in this material. However, our results 
indicate that the luminescence is a consequence of the large binding energy of the lowest exciton 
which stabilizes it against thermal relaxation. 

PACS numbers: 71.20.Nr, 71.35.-y, 73.22.-f, 78.20.Bli, 78.60.Lc 



Nanostructured forms of the semi-conductor M0S2 
have recieved much attention due to their potential as 
catalysts for desulferization of crude oil and more re- 
cently for (photo)-electrochemical hydrogen evolution [Ih 
[sl]. Bulk M0S2 is composed of two-dimensional sheets 
held together by weak van der Waals forces and indi- 
vidual sheets can be isolated by exfoliation techniques 
similar to those used to produce graphene(j]. Single lay- 
ers of M0S2 therefore comprise highly interesting two- 
dimensional systems with a finite band gap and have re- 
cently been proposed for nano-electronics applications [5|]. 

The optical properties of M0S2 have been thoroughly 
studied experimentally [6l4llj. The absorption spectrum 
shows two distinct low energy peaks at 1.88 eV and 



2.06 eV, which are denoted by A and B respectively [12 j. 
These features have traditionally been proposed to origi- 
nate from direct transitions between a split valence band 
and the conduction band at the K point of the Brillouin 
zone. Their Rydberg satellites, Zeeman splitting, and 
dependence on crystal thickness have been investigated 
in detail J9j|. Very recently, the quantum yield of lumi- 
nescence from M0S2 was shown to increase dramatically 
when the sample thickness was changed from a few lay- 
ers to a monolayer ^3J, [l^ and the results have been 
taken as an indication of a direct band gap in the single 
layer. This conjecture is supported by density functional 
theory (DFT) calculations, which indeed show a transi- 
tion from indirect to direc t gap w hen going from a few 
layers to monolayer M0S2 [iJ-llSJ. However, while DFT 
has proven very successful in describing the ground state 
properties of a large variety of materials, it is not ex- 
pected to provide a faithful description of excited state 
properties such as optical absorption and luminescence. 

In this letter we present first-principles many-body 
perturbation theory calculations for the electronic struc- 
ture of M0S2. We show that monolayer M0S2 has an 
indirect gap at the K — T transition, which is ^ 0.5 eV 
smaller than the direct gap at K . The experimental ab- 
sorption spectrum is reproduced by our calculations and 



it is shown that excitonic effects are very dominant due 
to the poor electronic screening in the two-dimensional 
layer. In particular, the A and B absorption peaks orig- 
inate from two bound excitons located at different posi- 
tions in the Brillouin zone and not from a split valence 
band originating from spin-orbit coupling and interlayer 
interactions as often stated in the literature p. This 
leads to a very different interpretation of the absorption 
spectra of M0S2 and related dichalcogenides [12| and re- 
solves the puzzle of why the double exciton is preserved 
in single layer structures, which do not have a split va- 
lence band 13|, [lj| . The high luminescence yield in mono- 
layer M0S2 can be attributed to the high binding en- 
ergy (0.7 eV) of the A exciton, which strongly reduces 
the probability that excited electrons decay to the low- 
est point in the conduction band before recombining. In 
bilayer and bulk M0S2 the difference between the direct 
and the indirect gap is larger and the exciton binding en- 
ergies are lower, thus enhancing the probability for non- 
radiative decay of the excitons. 

The many-body calculations described in this work 
were performed with the Yambo code[l7| and the ini- 
tial Kohn-Sham electronic structure was obtained with 



the DFT code ABINIT [l8l using the LDA exchange- 
correlation functional, a plane wave cutoff of 500 eV , 
and a 45 X 45 k-point grid centered at the F point. We 
have used the experimental lattice constants for M0S2 
and for the single layer calculations the repeated layers 
were separated by 23 A, whereas we used a separation 
of 29.3 A for the two layer calculations. For the one and 
two-layer calculations we used a Coulomb cutoff [19] in 
the direction perpendicular to the slabs to ensure con- 
vergence with respect to super cell size. The band struc- 
ture was obtained within the GqWo approximation (20| 
with a plane wave cutoff of 60 eV for the dielectric func- 
tion, which ensures convergence of the gap to within 60 
meV[28|. Unless stated otherwise the full frequency de- 
pendence of the response function has been taken into 
account. While the use of a plasmon-polc approxima- 
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FIG. 1: (Color online). Band structure of single layer M0S2. 
The Kohn-Sham band structure (blue lines) is seen to have a 
direct gap at K, whereas the QP band structure (red dots) 
has an indirect gap at the K — T transition. 

tion [2l| only induces a slight increase of the band gap 
of - 0.13 eV at the K point and 0.065 eV at the T 
point, these differences have consequences for the fine 
structure in the absorption spectrum. By assuming a 
static dielectric screening e(q, w) ^ £(q, 0), the Bethe- 
Salpeter Equation 22| can be recast into a two-particle 
excitonic Hamiltonian [23|, which is diagonalized on a 
basis of electron-hole pairs. In this way the excitonic 
eigenstates can be expressed as a linear combination of 
single-particle transitions 



|A)=5I^-fcl"c^)' 



(1) 



where v, c, and k denote valence band, conduction band 
and Brillouin zone wave vector, respectively. The absorp- 
tion spectrum is proportional to the imaginary part of 
the macroscopic dielectric function, which in the Tamm- 
Dancoff approximation can be written 



£2(0;) =27r lim v(q) \^ S(uj — E\) 

q-S-O ^ — ' 
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where E\ are the eigenvalues associated with |A). For 
a single valence and conduction band, the coefficients 
Ayck = A^k) can be regarded as an excitonic k-space 
" wavefunction" . 

In Fig. [H we show the quasi-particle (QP) band struc- 
ture of single layer M0S2 obtained within the GqWo ap- 
proximation as well as the LDA band structure. It is 
clear, that the inclusion of exchange and correlation at 
the Go Wo level, results in a severe modification of the 
band structure. The most prominent qualitative differ- 
ence is the minimum gap. Whereas the LDA band struc- 
ture shows a direct gap at the i^-point of 1.8 eV, the 



Go Wo band structure has an indirect gap between K and 
T of 2.3 eV, which is 0.5 eV smaller than the minimum 
direct gap at K. Very often, the true QP band structure 
can be well approximated by the LDA band structure 
subject to a "scissors operator" which shifts the conduc- 
tion bands by a fc-independent constant to reproduce the 
correct gap. While such an approach has proven useful 
for a range of different semiconductors, it would certainly 
be a poor approximation in the present case. 

The large direct band gap of 2.8 eV predicted by the 
Go Wo calculations may seem to be in conflict with the 
measured absorption spectrum of monolayer M0S2 show- 
ing peaks at 1.9 eV and 2.1 eV. However, the QP particle 
band structure does not include the Coulomb attraction 
between the electron and the hole involved in optical ex- 
citations. In materials with low screening, the electron 
and hole may form strongly bound states leading to exci- 
tonic features in the absorption spectrum well below the 
QP band gap. 

To include the effect of electron-hole interactions, we 
have solved the BSE and obtained the absorption spec- 
trum resulting from light polarized in the M0S2 plane. 
We have included the two highest valence bands and 
the four lowest conduction bands in the calculation (we 
have verified that the result does not change significantly 
when including more bands). The result is shown in Fig. 
[21 The inset shows the absorption spectrum calculated 
without electron-hole interactions. As expected from the 
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FIG. 2: (Color online). Absorption spectrum of single layer 
M0S2 reproduced from the experiments in Ref. [ij] (upper 
curve) and calculated from the BSE equation (lower curve). 
The calculated spectrum has been shifted down by 0.3 eV. 
The two peaks at 1.88 and 2.03 eV are identified with the 
A and B excitons. The inset shows the same BSE spectrum 
on a broader energy range along with absorption spectrum 
obtained directly from the Go Wo eigenvalues and LDA wave 
functions and with no electron-hole interaction included (red 
dashed curve). 
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FIG. 3: (Color online). Contour plot of the direct quasi- 
particle gap of monolayer M0S2. The gap shows two distinct 
minima which support two excitonic states. One exciton (A) 
is centered at the global minimum at the A'-point (QP gap 
~ 2.8 eV ) and the other exciton (B) is associated with a ring 
circling the F-point (QP gap ~ 3.1 eV ). The two minima are 
well separated by a gap barrier of ~ 0.6 eV. The contours 
are separated by ~ 0.1 eV and the gray dashed line is the 
Brillouin zone boundary. 



band structure it shows an absorption edge at '^ 2.8 eV 
corresponding to the direct QP band gap. However, 
when electron-hole interactions are included the absorp- 
tion edge is lowered to ~ 2.1 eV, which is well below the 
direct gap. The energy difference of 0.7 eV is the result of 
strongly bound electron-hole pairs. Usually, such strong 
excitonic effects are only found in insulators with low 
dielectric constants, but in the present case it is the two- 
dimensional nature of a single M0S2 layer which confines 
the electron-hole pair and gives the large binding energy. 
Such large excitonic effects have also been observed in 
hexagonal BN 2J, [25 1 and SiC [26J , which show similar 
two-dimensional confinement. The main panel in Fig. [2] 
shows the low energy part of the calculated absorption 
spectrum which has been shifted downward by 0.3 eV 
for comparison with experiment (black curve). It is seen 
to have four distinct peaks at low energies. The first of 
these is identified with the A exciton and peaks 2 and 
3 comprise the B exciton. In contrast to peaks 1-3, the 
fourth peak changes its weight and position when we in- 
crease the k-point sampling and although the peak seems 
to become less significant at high k-point samplings we 
have not been able to fully converge this feature in the 
spectrum. 

To demonstrate the nature of the peaks we can plot 
the magnitude of the expansion coefficients |>l^'^(/c)p 
appearing in Eq. [1] as a function of Brillouin zone co- 
ordinates. This is shown in Fig. [3] for peak 1 and peak 
3 (A and B excitons) along with contour lines indicating 
the direct QP gap. Peak 2 has nearly identical weights as 
peak 3 but has an additional small weight near K . The 



QP gap has a global minimum situated at the K-point 
and a well defined local minimum located in a ring de- 
fined roughly by |A;| ^ \K\/3. These minima are seen to 
support the A and B excitons. The expansion coefficients 
also allow us to calculate the exact binding energy of the 
two excitons. Using that the exciton energy without the 
e-h interaction is E'^% = J2k l^^'''(fc)l'(e?/ " ^vk) ^e 
find that the binding energy of the A and B excitons are 
0.7 eV and 0.8 eV, respectively. The slightly larger bind- 
ing energy of the B exciton agrees with the fact that it is 
more delocalized in the Brillouin zone and thereby more 
localized in real space resulting in stronger coulomb in- 
teraction between the electron and hole. Equivalently, 
we may say that the A exciton has a smaller hydrogenic 
binding energy due to its smaller effective mass result- 
ing from the high "gap curvature" at the K-point. In 
fact, we may estimate the value of the exciton binding 
energy as E^^ — m^/2e1 where E^^ is the binding en- 
ergy in atomic units, ?7ie the exciton effective mass, and 
£r the dielectric constant. We have calculated £r within 
the random phase approximation and obtain e|| = 4.1 
and £_L = 1.4 for parallel and orthogonal polarization 
respectively. The exciton effective mass at K is 0.42 
which gives E^^ — 0.32 eV for the parallel component 
of er- However, since the exciton is confined in a nearly 
two-dimensional region, it may be more reasonable to 
consider a two-dimensional hydrogen model, which gives 
El^ = AEl^ = 1.3 eV [13]. While, such estimates are 
not reliable for detailed quantitative analysis, the hydro- 
gen model does give the right order of magnitude for 
the binding energy and the BSE value of 0.7 eV is close 
to the average value of the two-dimensional and three- 
dimensional models. 

The recently measured high luminescence yield in 
monolayer M0S2 13|, |lJ|, seem to be in conflict with 
the Go Wo band structure shown in Fig. [TJ This is be- 
cause the indirect gap suggests that electrons excited at 
K would relax to the T point by phonon scattering from 
which point they would not be able to recombine directly 
with the hole left at K. This picture is, however, based 
on the QP bands and do not take the strong binding en- 
ergy of the electron- hole pairs into account. Indeed, the 
energy of the A exciton is 0.2 eV lower than that of a free 
hole at K and a free electron at T (since the indirect gap 
is only 0.5 eV below the direct gap). Even if the bind- 
ing energy of the indirect exciton, consisting oi a K — T 
electron- hole pair, is comparable to that of the direct 
exciton at K, luminescence could still take place if the 
nonradiative relaxation of the exciton is slow compared to 
the recombination time. This is not an unreasonable as- 
sumption since the downhill relaxation of the exciton re- 
quires that the coherence between the electron and hole is 
conserved during the electron-phonon scattering events. 
Therefore, we propose that the high luminescence yield 
in monolayer M0S2 is the result of strong excitonic bind- 
ing energies and associated low relaxation rates, rather 




of many-body perturbation theory calculations. 



FIG. 4: (Color online). Left: Band structure bilayer M0S2. 
Right: Band structure of bulk M0S2. The blue lines are 
the Kohn-Sham band structure and the dots are Go Wo band 
structure. 



than a direct band gap as previously suggested 13lll4|. 

For comparison we have also calculated the band struc- 
ture and absorption spectrum of bulk and bilayer M0S2. 
We find quantitative differences from monolayer M0S2 
which are expected to lower the quantum yield of lumi- 
nescence. To illustrate this, we show the band structure 
of bulk and bilayer M0S2 in Fig. 2] The results were ob- 
tained with a 30 X 30 k-point sampling, and the bulk band 
structure was calculated using the plasmon pole approx- 
imation. For the bilayer (bulk), the indirect gap is 0.8 
eV (0.9 eV) lower than the direct gap. The same value 
is only 0.5 eV for the monolayer. Furthermore, in the 
bulk material the screening is much larger than in mono- 
and bilayer structures and the exciton binding energies 
are less than 0.1 eV. Consequently, the phonon mediated 
scattering of the photo excited excitons from the direct 
to the indirect gap should be more effective in bilayer 
and bulk structures in agreement with the fact that no 
or very low luminescence is observed in these materials. 

In conclusion, we have performed first-principles cal- 
culations on M0S2 and shown that both monolayer and 
multilayer structures have an indirect band gap. We pro- 
pose that the recently observed high quantum yield of 
luminescence in monolayer M0S2 can be attributed to 
the large exciton binding energies leading to low thermal 
relaxation rates. Our results show that the quasiparti- 
cle band gap obtained with the GqWo approximation is 
qualitatively different from the LDA band structure and 
leads to a very different interpretation of absorption and 
luminescence spectra. In particular, the A and B peaks 
in the absorption spectrum are not associated with a va- 
lence band split by interlayer interactions but rather is 
due to two distinct, strongly bound excitons. 
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